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ABSTRACT 

The paper is concerned with a theoretical study of shearing flow 
bounded by a wavy surface with consideration of turbulent flow above 
the air-sea interface. Account is to be taken of turbulence in the flow 
through the use of Reynolds stresses associated with turbulent flow. 

An adaptation of the mixing length theory as applied to pipe flow is made 
for channel flow and the resulting mixing length versus energy relation- 
ship is incorporated in the Reynolds stress term. Finally, the rate of 
wave growth is calculated from the normal surface pressure in phase 
with the wave slope. 

A curvilinear coordinate system which follows with the wave train 
is used in order to simplify the formulation of the problem. All para- 
meters are non-dimensionalized and the analysis is made considering 
a velocity profile adapted from pipe to channel flow. 
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non-dimensional wave celerity 
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1. INTRODUCTION 



A semi-empirical study of turbulent air flow over a wavy surface 
is made in which normal and tangential stresses on the surface are 
examined with the object of estimating the growth of water waves due 
to a turbulent wind . 

Ursell (1956) gives a summary and a critical review of the theories 
of wind-wave formation extant in 1956. The review gave the necessary 
impetus for further investigations, both theoretical and experimental, 
in this area. Since 1956 there has been much work done on the generation 
of water waves by surface stresses induced by the wind. Miles (1957) 
and Benjamin (195 9) considered laminar air flow over a wavy water 
surface and deduced formulae for exponential wave growth due to an 
instability mechanism. Later Phillips (1966) included air-stream 
turbulence in the instability model from a semi-empirical point of view. 

In his initial work, Miles (1957) assumes an inviscid, incompress- 
ible, two-dimensional, parallel shear flow coupled with a prescribed 
two-dimensional deep-water gravity wave. Turbulent fluctuations in 
the perturbation equations of motion are neglected, even though they 
are decisive in maintaining the mean shear flow. Furthermore, the wave- 
induced velocities and pressure are small, which justifies linearization 
of the equations. The water is also assumed to be inviscid, incompress- 
ible, and irrotational , and the slope of the displaced surface is assumed 
to be small. Only that component of the aerodynamic force in phase 
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with the wave slope is important in the wave-generation mechanism. 

The viscous drag forces of the air are assumed to be negligible compared 
with the normal pressure in their effect on the surface wave. 

This model is used to determine a first-order approximation to the 
disturbed motion of the air and the consequent energy transfer to the 
wave. The velocity distribution is taken to be that for turbulent flow 
over an aerodynamically smooth water surface. The result rests on the 
solution of the inviscid Orr-Sommerfeld equation which has a singularity 
at the height where XT = c (the critical level) . An exponential 
increase in amplitude of water waves with phase velocity c occurs if 
the mean air flow velocity U (y) in the direction of wave travel varies 
with height above the surface such that TJ at the critical level is 
negative. This is almost always the case. 

Phillips (1966) attempts to include wave-induced turbulence in the 
instability model by including in the momentum flux equation an 
appropriate integral over the entire layer of air of the contribution due 
to perturbation eddy viscosities. This integral is multiplied by an 
unknown constant A which is estimated from experimental results from 
flow over fixed wavy surfaces. The uncertainty in A is estimated to 
be + 50 per cent. Phillips concludes that even though there exists 
a high degree of uncertainty in A the results do provide a framework 
in which experimental measurements may be interpreted. 

According to Miles (1967) , experimental results indicate that a 
theoretical model based on laminar flow may be adequate in the 
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laboratory but not on an oceanographic scale. The inviscid, laminar model 
underestimates the energy transfer from wind to waves over at least a 
significant portion of the spectrum for an open sea. This suggests that 
the significance of wave -induced perturbations in the turbulent Reynolds 
stresses for momentum and energy transfer from wind to wave must 
increase with an appropriate scale factor. Miles concludes that further 
theoretical work in this area appears to demand some 'ad hoc 1 hypothesis 
to describe these perturbation Reynolds stresses. 

This can be done by considering a slightly generalized form of the 
classical mixing-length hypothesis. A means must be found to predict 
the value for the change in mixing length upon which the perturbation 
Reynolds stresses are based. Zagustin (1963) has given a semi-empirical 
solution for turbulent flow in a pipe which gives good agreement with 
experimental data and which predicts the value of the mixing length. 

This model is generalized below to include flow over a wavy surface in 
a manner analogous to that of Miles and Benjamin in that it assumes 
the flow to follow the wavy surface. Viscous stress terms are replaced 
by the Reynolds stresses associated with turbulent flow. Normal Reynolds 
stresses are neglected. Zagustin' s method of closing the system of mean 
turbulent equations of motion is not unique (see, for example, Squire, 

1959) , but there seems to be no good reason to use any other model in 
its place. The Investigation below details the application of Zagustin' s 
work to this model and outlines the solution of the differential equations 
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using numerical methods and appropriate boundary conditions. The 
unknown parameters in the model are fitted to experimental data describing 
the normal pressure in phase with the wave slope. 
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2. FORMULATION OF THE MODEL 



The zero-order flow is taken as flow in a two-dimensional channel. 

All lengths are made non-dimensional using the channel halfwidth, H. 

Velocities are made non-dimensional using the midstream velocity, TX • 

Z 



All stresses are made non-dimensional by dividing by o \J 0 

Ja~ 

is the density of air. 



, where 



2 . 1 THE ZERO-ORDER FLOW 



The velocity is expressed in terms of a stationary and a fluctuating 
(turbulent) part: 



( 2 . 1 ) 



iu. 



= IT 



4- WL 



Following Zagustin, the Reynolds shear stresses ll are written as 



( 2 . 2 ) 



Y - L (cjro-cLc^) 



where 0^ is the mean speed. This can be considered a defining 
equation for L_ (although here it appears as the usual mixing length) 
in that L will not be prescribed, but will be part of the solution. Then 
the two-dimensional mean equation of motion is 

(2.3) ( ¥ • grad) U = -grad P + |K x grad 7T 

A 

where IK is a unit vector normal to the plane of the flow, P is the 

I ^ 

pressure, and K = |_ grad q is the eddy viscosity. The turbulent 
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\ I 

kinetic energy, E' = <iL* IU. , is assumed to have the form 
(2.4) E ' = C L 2 (grad q) 2 



where C is an unknown constant. It is assumed that the rate of 
dissipation of E' is proportional to E'. 

A flux of turbulent kinetic energy into any volume element is assumed 
to occur when the quantity I— varies over the element. Thus, the 
rate of inflow of E' is 




where n is the normal to any surface and is an unknown constant 

depending on the turbulent intensity. The turbulent kinetic energy is 
constant at any point, and the above assumptions then imply that 



( 2 . 6 ) 



(grad 



L 



n)dS 



+ 




E' dV = 0, 



where B is an unknown constant. Applying the divergence theorem, 
considering V to be an arbitrary volume, and using equation (2.4), 
equation (2.6) becomes 



(2.7) 



V 2 L + B L 2 (grad q) 2 = 0 . 



Equations (2.3) and (2.7), together with the continuity equation, 
form a closed set for the mean velocities and the quantity L_ . The 
constant B must be determined experimentally. 
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For flow between two parallel walls as shown in Figure 1, the 



continuity equation is trivial, and equations (2.3) and (2.7) become 



( 2 . 8 ) 



-grad P + 



d_ r l? 

dy L 



V 



dy 



= 0 



and 



(2.9) 



2 L + B [_ 2 (d ~P~ 



dy^ 



dy 



= 0 . 




y = 1 



y = -1 



Defining B = A for the zero-order flow case and setting 
2 



= | grad P | , 



equation (2.9) becomes 

,2 



( 2 . 10 ) 



d‘ L , A 

dy 2 2 



I 2 (d TT 

L dy 



= 0 , 



and equation (2.9), after integrating, becomes 



(2.11) 



hy = |_ 2 teL 



V 



dy 



Equations (2.10) and (2.11) give 



( 2 . 12 ) 




19 



which has the boundary conditions: 



L 



dy 



= 0 and 



L - Lo. 



L 



y = i 



y = 0 

where I o is a quantity analogous to the roughness length in mixing- 

length theory. Integrating equation (2. 12) twice and applying these 
boundary conditions results in 



(2.13) 



L - itu -v 3 > + L, 



u 



Next, solving equation (2.11) for d V and substituting for 

dy 



L 



gives 



(2.14) 



XT 



dy 



L 






+ ^-(1 - y 3 ) 

o 1Z 



By defining 



b = 1 + 



12 



.L 

Ah 



and 



z = 



. _z_ 



V7 



equation (2.14) simplifies to 



(2.15) 



12 

A 'J hb' 




y~z~* dz. 
1 - z^ 



3 

Integrating equation (2. 15) , applying the no-slip condition at z = - \f 
and returning to the y coordinate gives 
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(2.16) 



TJm - 



8 



A 7 hb' 



£ tanh * GU - tanh 



TiT 



7F 



The mean velocity at mid-channel is 



T7 (y) 



= l 



y = 0 



which, when applied to equation (2.16), results in the expression 



(2.17) 



8 



A VET? 



tanh 



-1 i _ 



(7F) 



= 1 



Equation (2.16) therefore becomes 

TJm - ■ ■"'" 1 ( ' 



tantivy 



(2.18) 



1 - 



/F 



tanh” 1 ^! 



L 



The expression (2.13) for l_ can be simplified by recalling that 

von Karman's constant, K , gives the rate of growth of I at the 

wall. Therefore, 

iL 



(2.19) 



K = 



dy 



y = - 1 

where VC ^ 0.4 is the accepted value. Then, from equations (2.13) 
it is found that 



( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 



Ah = 1.6 



b =1 + 7.5 



L 



L - 



b - V 3 

7.5 
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Following Benjamin (1959), the equations of motion in terms of j| 



and are 



(2.30a) 

b) 



KOi, C- ) + i/2J^(0| 2 + Ci, 2 ) 

i ( 'H^ + ’ + 1/2 V^ 2 + K ? 2) = -P> 1 + Z1 



where the subscripts denote partial differentiation. 

The stream function is composed of a zero-order term plus a 
periodic perturbation. When there is no wave, the stream function is 



(2.31) ^ 




where + y = 1 . 

(J ( Vp is the velocity profile from equation (2. 18) expressed in 
curvilinear coordinates 



(2.32) 



TP* 



tanh 



)= 1 - 



-1 [y? (1 ~ t?)3/2 1 . 

«nh-l (i.) 



When the water surface is disturbed by a stationary wave, the stream 
function becomes 



(2.33) ^ ,>J) = (^ D 





where F(>| ) is to be determined. 

The mixing length is expressed in a form similar to that for : 



(2.34) 



Ji = 
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L 



X 



where L (T, ) is given above and 4 ) is to be determined , 

Similarly, 



(2,35) 



B = 4 + a c lk ) 



where A is given above and ^ is an unknown parameter to be 
determined empirically. The velocity components parallel to ^ and 
are given by 

(2.36a) u=J 1/ ^= U - c + a |V(?{) + U (>j )e" kl ?j- e ik ) 



b) v 1 



1 '2 






= -ika 



-T F(>j) + -o] e ik l 



To formulate the differential equation in to be solved, 

substitution is made for , B and q in equation (2.7) expressed in 
curvilinear coordinates: 



(2.37) 






+ Bl ^-) 2 = 0 . 



Simplification by elimination of all terms of 0(<- 2 ) and above yields 

(2.38) X ♦ i. (AE k 2 ) + AHF" + ^ D + 

Ae" k> ? (G - kD) = 0 



where 



D = IT 2 u 
B - IT 2 L 



and 
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Similarly, when equations (2.30) are expanded, cross differentiated 
to eliminate pressure terms, and linearized in a, one obtains 



(2.39) 



iL 

2 



^ ( \J - C)(F" - k 2 F) - X}" F sgn( U' ) = 



+ k 2 )(HF" + E-f ( ) + 

d>| 

e~ kY { (G" - 2kG' + 2k 2 G) 



where 

g -U'U"L 2 

H = U'L 2 

and where "U” is positive for the lower half of the channel. Note that 
a singularity occurs at mid-channel ( ' = 0) . Equations (2.38) and 

(2.39) are to be solved for F and . 

The boundary conditions are as follows. At the interface where 
V = o, XJ (0) = 0 to satisfy the no-slip condition, and v = 0 because 
the wave is stationary. Therefore, from equation (2.36b): 

(2.40) F (0) = c. 

The velocities in the water must take on negative values in order that 
the shearing stress across the interface be continuous. The variable 
component of the tangential water velocity associated with an irrotational 
wave is of the form 



(2.41) 



u' = a 
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The total tangential water velocity is 



(2.42) u = -c - a yr / <&_ , n e ik £ . 



IL 



Defining the Froude number: 

TT 2 

(2.43) H" = % T ° , and invoking the no-slip condition, the 

11 gH 



tangential air velocity is 
(2.44) 



(2.45) 



u = -c - a 



u = -c - a 



y^p e ik > at =0, or 



0 ^ 



at = 0 



where 



p-IF 



. From equation (2.36a) the second lower boundary 



condition then becomes 



(2.46) 



F'(0) = - TT'(O) + j?> . 

The perturbation in Si at the water surface, Jl t (0), must also be 
fixed. This must be regarded (along with ) as an unknown parameter 
in the problem, to be adjusted empirically. 

The perturbed flow must die away at some large distance (relative 
to k ^) from the wavy surface. Thus, as VJ becomes large, 

— t- (c -TJ )e~ k1 1 

- U' e" kl I 

0. 



(2.47) 


F(V|) 


(2.48) 


F’(7j ) 


(2.49) 
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2 . 3 NORMAL PRESSURE , SHEAR STRESS AND WAVE GROWTH 

The pressure on the wavy surface is the sum of that due to normal 
atmospheric turbulence and that generated by air moving over the 
undulating surface. The latter pressure predominates in the generation 
mechanism considered here and is expressed as 

(2.50) P(0) - a( %/• + icf )e ik 

where % and & are to be determined . Differentiating equation 

(2.50) with respect to l» , substituting the result into equation (2.30a), 
and simplifying yields 

(2.51) P(0) = [TT’f - ( U " c)F’ - (HF" 

+ E^,) + (G‘ - Gk) 1 

J y\ = o 

which gives '’X' and & . 

In the same manner, the shear stress at the wall is expressed as 

(2.52) (0) = a( 6 + i E )e* k ^ 

where £ and £, are to be determined. To evaluate the shear 
stress, equation (2.2) is expanded in terms of ^ : 

( 2 . 53 ) T = JLh J 2 

Substituting for JL , j, (+/ , or their derivatives where required 

obtains 
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(2.54) 



T (0) - ( 2 LU'[LU' + L< f " + T7 M) | 



which gives 6. and £. . The shear stress seems to be relatively 

unimportant to wave growth, however, and will not be considered herein. 

To compute the wave growth rate the normal pressure in phase with 
the wave slope is required. Expansion of equation (2.50) and comparison 
with equation (2.25) for the wave form show that a positive <£ results 
in wave growth. The average energy input due to normal pressure is 



(2.55) 








'o 

or, in dimensional terms. 



d 



* 



t a 2 TT 2 r 

(2.56) E = 2H W ?o- Uo 

o 

where CO = gk is the wave frequency. 
The surface energy of the wave is 



(2.57) 



E = 1/2 



fur 



ga 



where 






is the water density, so that 



(2.58) 



£ = 



da 



fa 9a d, 



Equating equations (2.56) and (2.58) and solving for gives 



(2.59) 



da - a- 
dt 






2gH 



to 
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Integrating produces the familiar exponential wave growth equation: 



( 2 . 60 ) 



a = a 0 exp 



LQ*_ 



IF £ up t 
2 



~\ 



>** - J 

The rate of energy input due to surface shear is 



( 2 . 61 ) 



E = h 



7 \ 






u d^ 



/ 



or 



( 2 . 62 ) 





This energy is expected to go, at least initially, into a surface shear 
flow. 
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3. SOLUTION OF THE FIRST-ORDER PROBLEM 



A solution for the function F and its derivatives must now be 
found so that values for normal pressure, shear stress and wave growth 
rate may be determined. There are several numerical integration methods 
available for this purpose, having the common property of starting at 
one boundary and ending at the other. However, the problem as 
formulated has conditions at both boundaries, and a method of super- 
position of solutions must be used. 

Equations (2.38) and (2.39) can be rewritten as six linear first- 
order ordinary differential equations. That is letting: 



(3.1a) F = Y 1 



b) F' = Y 2 



c) F" =Y 3 



d) F" = Y 4 



e) 




f) 




gives six ordinary differential equations: 



(3 . 2a) Y' = Y 2 



b) r - y 3 

0 r - y 4 
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d) VH{f[ ( f - <=) (Y 3 - *\) - n" Yl j - 2H'Y 4 

- (H M + k 2 H)Y„ - eJL! - 2E’Y c - (E" + k 2 E)Y c 

o bo 

- e~ k "^ (2Gk 2 - 2G'k + G")j* 

and 

e) Y’ = Y 

5 6 

f) Y' = -Y c (AE - k 2 ) - AHYo - )T D 

6 5 J 

-Ae" k ^ (G - kD) . 

There are now six equations, each requiring an initial input at the lower 
channel boundary. Linearly independent solutions are found by starting 
repeatedly at the lower boundary with different initial values, Y(0) . The 
final solution is obtained by superposing the independent solutions to 
fit the required conditions at the upper boundary. 

The independent solutions, denoted by the superscript j, have 
initial values given by 

(3.3) Y^ +1 (0) = ^ j at 1^1 = 0, 

where is the Kronecker delta. The solution to the homogeneous 

form of equations (3,2) can be written as a sum of independent solutions 

to the homogeneous equation: 

6 

(3.4) Y. = ^ J, Y} 

1 J 1 
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where the J's are complex constants . Since 
F (0) = F'(0) = Jl % (0) = 0 

for the homogeneous system, the constants Jj, J 2 , and J,. in equation 
(3.4) must equal zero. The complete solution is achieved by adding to 
the homogeneous solution a particular solution of equation (3.2). For 
this particular solution, the initial conditions are 

(3.5a) Yj = c 

b) Y p = £ - U 

c) Y p = 0 

3 

d) Y p = 0 

4 

e) Y p = Ji t (0) 

f) Y p = 0. 

6 

Then the complete solution, satisfying all the boundary conditions at 
^ = 0 , is 

(3.6) YjOj) = J 3 y3(>P + J 4 Y 4 (>() + J ) + Y^(^), 

where the J's are to be determined. 

The position of the upper boundary, where conditions (2.47), (2.48), 

and (2.49) hold, is not known, and must be determined by a convergence 

criterion. The solution is determined by applying the upper boundary 

condition at various increasing values of = d, solving the resulting 

equations (3.6) evaluated at d for the J's, and evaluating the wall 

pressure each time using equation (2.51) . When the wall pressure 
determined in this manner converges , the procedure is stopped. 
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4. THE COMPUTER SOLUTION 



The solution for the function F and its derivatives has been 
obtained by solving the problem as outlined in section 3 on an IBM 360- 
67 computer using complex double-precision arithmetic. Major sub- 
routines are incorporated in the program to calculate the velocity profile, 
perform the integration using the Runge-Kutta Gill fourth-order method, 
compute the shear stress, calculate the wave growth rate and evaluate 
the stream function. Appendix I contains the format of the program and 
a set of results obtained. 

The accuracy of the Runge-Kutta Gill integration method depends 
on the size of the step used in the iteration. The value 10 was 
selected in order to determine whether the surface pressure converged 
on a particular value. Because surface pressure did not converge for 
this step size, it was decreased until convergence was achieved. A 
step size of 10 was necessary to obtain convergence of the wall 
pressure, and the program running time was about 112 minutes. A final 
optimum step size of 10 ^ for less than 0.3 and 2.5 X 10 ^ for ^ 

greater than 0.3 was then adopted. This step size gives convergence 

„ r 

on a result essentially the same as that for 10 and a running time 
of 65 minutes. 
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5 . RESULTS AND CONCLUSIONS 



The zero-order velocity profile is given by equation (2.18) which 

has 1 (0) as a parameter. A plot of U (y) for several values of l — (0) 

appears in Figure 3 . From the figure it can be seen that as the parameter 

| (0) decreases, the shape of the profile changes from nearly linear to 

almost a step function. A representative logarithmic profile is also 
shown in the same figure so that the two profile shapes may be easily 
compared . 




Figure 3 . 
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The remainder of the analysis is concerned with the determination 
of suitable values for the parameters |_ (0) , ^2, (0) , and r . For 

the channel under consideration: 

H = 10 meters 
U- 10 meters/second 

which results in a Froude number of 1.018. For the wave number, 
k = 15, 

the instability mechanism considered is important. These values for 
the Froude number and the wave number have been used throughout the 
analysis . 

The following parameter values were first used: 

| (0) = 2.22 X 10~ 2 

JE, (0) = 0 

X' = 10 . 

The wall pressure failed to converge for these values. Convergence 
was assumed to have resulted if wall pressure values were nearly 
constant over a significant range of d (say 0.1) or reached a minimum 
value at some point in the channel. (This latter criterion is not 
completely defensible, but the singularity in equation (3. 2d) at = 1, 
itself a peculiarity of the Zagustin assumptions, is probably responsible 
for the erratic behavior of the wall pressure as was then 

varied from -1000 through zero to +1000 and given imaginary values as 
well. The wall pressure changed only in the fourth or fifth decimal 
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which was considered insignificant, although changed markedly. 

Thus, the solution seems insensitive to ^ . 



New values for the parameters were then used: 

L (0) = 5 X 10~ 5 

Jl, ( 0 ) = 0 
& = 10 . 

The wall pressure converged for these values, but the value was too 
high in magnitude and of the wrong sign. A value nearer to 5 (Miles, 
1957) was considered appropriate. 

With 

L (0) = 5 X 10" 5 

JP, (0) = 1 

=10 



the resulting wall pressure was positive but still too high in magnitude. 



Figure 4 summarizes the results of varying 





( 0 ). 




Figure 4 . 
3 7 



It appeared that the wall pressure was more directly dependent on 
(0) than on any other parameter. Subsequent investigation showed 
that for any >| the wall pressure was a linear function of the real or 
imaginary part of (0) providing that the other component was held 
constant. For 

I (0) = 5 X 10” 5 

i/O, = 0.65 788 

= 10 

the resulting minimum wall pressure was 3.83 at d = 0.48. Figure 5 
shows the resulting wall pressure in this case versus d. 

$ 

1.0 

0.8 

0.6 

0.4 

0.2 




Wall Pressure, 
Figure 5 . 



was then given a value: 

ft* = 1000 . 
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No significant change occurred in the wall pressure. Convergence, as 



defined above, always occurred near d = 0.5. 



Table 1 summarizes the values of the wall pressure, shear stress. 



and 


wave growth rate versus , 


(0) at d = 0. 


5 for parameter values: 








L ( 0 ) = 5 X 10 5 










r = io . 






2 






Out-of-phase 


In-phase 


Wave Growth 


(0) 




Wall Pressure ,<£ 


Shear Stress 


, £ Rate , .Bl. X . 












o 0 O C 

)o3 2~ 


0.0 


+ 


Oi 


-6495.3 


-29.78 


-1.61 


0.2 


+ 


Oi 


-4519.4 


-14.89 


- 1.12 


0.4 


+ 


Oi 


-2543.6 


0.001 


-0.63 


0.6 


+ 


Oi 


-567.7 


14.87 


-0.141 


0.8 


+ 


Oi 


1408.2 


29.76 


0.349 


1.0 


+ 


Oi 


3384.0 


44.64 


0.838 


0. 65788+0i 


3.83 


19.19 


0.001 


0.0 


+ 


0.21 


-6495.3 


-29. 78 


-1.61 


0.0 


+ 


0 . 4i 


-6401.8 


-29.78 


-1.58 


0.0 


+ 


0. 6i 


-6308.3 


-29.78 


-1.55 


0.0 


+ 


0 . 8i 


-6214.8 


-29.78 


-1.52 


0.0 


+ 


l.Oi 


-6221.3 


-29.78 


-1.49 


0 . 65788t 0 . li 


50.9 


-29.78 


0.013 




Table 1. 


Wall pressure, 


shear stress, 


and wave growth rate 








as a function of the perturbation surface roughness. 



39 



It can be seen that imaginary values for (0) have no effect on the 
shear stress. The required value for (0) is 0.65788 + Oi for a wall 

pressure of 3.83, but the shear stress seems high. (However, the shear 
stress value is a first-order contribution to the total shear stress and is 
reduced by the fractional value of the wave amplitude.) Exponential 
wave growth rates, due to the wall pressure alone, appear reasonable. 

To reduce the shear stress to a more acceptable value in relation to the 
normal pressure, parameter values 

L(0) = 5 X 10" 5 

J?, (0) = 0.4 + 5.45i 
^ = 10 

were used, and give at = 0.5 

Normal pressure = 5.0 
Shear stress = 0.001 
Growth rate = 0.001. 

However, the wall pressure failed to converge, and had the form shown 
in Figure 6. 

Equation (2.33) for the stream function has been evaluated with 
the values of F calculated from the constants J determined at d = 0.48 
where the wall pressure was found to converge. Parameter values 
were 

L(0) = 5 X 10~ 5 

i.(0) = 0.65788 

= 10 . 
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Wall Pressure 



Figure 7 shows a plot of the streamlines with the 'cats eye' present 
and centered about 30 degrees ahead of the wave crest. The 'cats eye' 
or closed loop is indicative of the wave induced perturbations in the 



flow and characterizes the region of turbulence (Phillips, 1966). Higher 
in the channel the streamlines become less wavy and approach straight 
lines. Table 2 shows values of the stream function near the critical 
layer. 

In view of the foregoing limited results, the following tentative 
conclusions can be made: 

1) The inverse hyperbolic tangent velocity profile, a result of 

the theory of the model, can be fitted to almost any known profile type. 
The profile which has a large value for -u near the critical level 
and for which the parameter |_ (0) has the value 5X10^ produces 
surface pressures which can be made to converge on the desired value 
by varying the wall perturbation in . 

2) Wall pressure is extremely sensitive to very small changes in 

-^i(O) and insensitive to moderate changes in ^ . 

3) Reasonable wall pressure (and wave growth rates) can be obtained 
with an appropriate choice of (0). This indicates that small, unequal 
variations in the surface roughness, such as would be caused by 
parasitic capillaries, are very important in the wave-generation process. 

4) Shear is independent of imaginary values of (0). Therefore, 



individually or in combination as indicated by further analyis of the model. 




X or 




must be varied 
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Streamline Pattern 



Based on the above results, the following suggestions are made: 

1) 1 (0) should be increased slightly to 

1 ( 0 ) = 10" 4 

in order to reduce the shear stress relative to the wall pressure. One 
computer run indicates that the shear stress is halved for this value. 

2) Additional wave numbers should be investigated, say, in the 
range 5 k C 20, and their effect on wall pressure and shear stress 
noted. 

3) The model should be checked using data from Stanton and 
Motzfeld as quoted in Ursel! (1956), and with recently obtained laboratory 
data (see, for example, Zagustin, et al . , 1966). 

In summary, the model appears to be fairly well suited to a 
theoretical study of air-stream turbulence and wave generation. 
Unfortunately, the long time required for a computer run prevented the 
calculation of further data using appropriate parametric values . 
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k l 





0.0 


0.1 


0.2 


0.3 


0.4 


7 

0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.002 


-1.522 


-1.256 


0.508 


0.434 


0.216 


0.005 


-0.683 


-0.571 


-0.227 


0.218 


0.5 94 


0.008 


-0.299 


-0.253 


-0.080 


0.154 


0.359 


0.010 


-0.138 


-0.116 


-0.007 


0.145 


0.285 


0.015 


0.110 


0.101 


0. 129 


0.183 


0.243 


0.020 


0.244 


0.229 


0.238 


0.269 


0.309 


0.025 


0.325 


0.315 


0.336 


0.381 


0.432 


0.035 


0.424 


0.440 


0.525 


0.646 


0.758 


0.045 


0.506 


0.556 


0.722 


0.940 


1.128 








k l 








0.5 


0.6 


0.7 


0.8 


0.9 


*7 












0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.002 


1.533 


1.266 


0.518 


0.426 


-1.522 


0.005 


0.757 


0.645 


0.300 


-0.145 


-0.520 


0.008 


0.456 


0.410 


0.236 


0.003 


-0.202 


0.010 


0.357 


0.336 


0.228 


0.075 


-0.065 


0.015 


0.284 


0.293 


0.265 


0.211 


0.152 


0.020 


0.343 


0.359 


0.349 


0.319 


0.278 


0.025 


0.470 


0.480 


0.459 


0.414 


0.363 


0.035 


0.817 


0.801 


0.716 


0.594 


0.482 


0.045 


1.212 


1.162 


0.996 


0.778 


0.591 



Table 2. Values of the stream function near the critical layer 
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DENOM = A2 - A3R 
U1 = (CON*A3) /DENOM 

INPUT VALUES AT ETA = 0. FOR I ^HOMOGENEOUS EQUATION 
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INHOMOGENEOUS EQUATIONS 

6 C 1 EL ( 2 ) = -Y( 5 )*(A*E - WK ? ) - A*H*Y(T) - GAMMA*!) 



A*CDEXP(-ETAK)*(G - WK*D) 

) = ( l./H)*(Z*WK/ 2 .*((U-C)*(Y( 3 ) - WK 2 *Y( 1 )) - tJ 2 *Y(l)) 
.*H 1 P*Y( 4 ) - ( H 2 P + WK 2 *H)*Y< 3 ) - E*EL( 2 > - 2 .*E 1 D *Y( 6 ) 

E 2 P + WK 2 *E)*Y( 5 ) - COE XP ( - ET AK ) * ( 2 . *G*WK 2 - 2 .*WK*G 1 P + G?P) 
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MODIFIED JORDAN ELIMINATION METHOD 
IMPLICIT REAL *8 (A-H, O-Z) 
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